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Abstract 

Although two pion decays of K, D and B have similar isospin structures, there are dramatic 
differences in the ratios Rk and Rd,b of amplitudes from AI = 3/2 and AI = 1/2 interactions. 
In K — > 7T7T decays there is the famous AI =1/2 rule with Rk ~ 1/22, whereas in B(D) — > 7T7t 
decays the ratios -Rd,b are of order one and therefore there is no such a rule. In this work we study 
decay amplitudes in B(D) — > 7T7t using QCD factorization calculations paying particular attention 
to discrepancies between data and theoretical estimates. Since isospin does not play a special role 
in factorization calculations, no AI =1/2 rule is expected. We find that theoretical results on 
the size of the amplitudes are in qualitative agreement with data. However the phases for the 
amplitudes are very different. We show that the effects of re-scattering between the two pions in 
the final state can play a crucial rule in understanding the differences between B(D) — * irir and 
K — > 7T7T decays. We also comment on the role of isospin analysis which applies to the study of CP 
violation in B — > tttt decays. 
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It is well known for many years that there is a AI = 1/2 rule for K — > 7r7r decays. 
In these decays the ratio Rk of the isospin 1 = 2 amplitude for the two final pions from 
AI = 3/2 interaction to the 1 = amplitude from A J = 1/2 interaction is about 1/22 
which is much smaller than naive expectations. Why there is such a big difference between 
the two amplitudes presents a serious challenge to theory. Many attempts have been made 
to explain the AI = 1/2 rule. Since the pioneer papers by Gaillard, Lee and Altarelli, 
MaianiQ, it was hoped that the AI = 1/2 rule may be explained by the running of the 
relevant Wilson coefficients. However, the different evolution of the Wilson coefficients of 
the operators which correspond to A J = 1/2 and AI = 3/2 respectively can only result in a 
factor of 2 ~ 4, as the data demands a deviation of more than 20, another factor of 6 ~ 10 is 
still missing. Other theoretical scenarios have been proposed to explain the difference. For 
example, the final state interaction Q, the contribution of the penguin-induced operator Oq 

and even the diquark intermediate states j^j]. 

Since the mesons K, D and B all have the same isospin structure, it is natural to expect 
that in two pion decays of B and D there are similar AI =1/2 rules just like in K —>■ titt 
decays. It is interesting to study in more details the isospin amplitudes of B(D) — > tttt 
decays. Moreover, one hopes that the study of whether the AI =1/2 rule holds or not in B 
and D cases can provide crucial information to understand the mechanism which enforces 
the AI = 1/2 rule in the K case. 

The decay amplitudes for K, D and B into two pions can be parameterized in the same 
way as, 



where P is one of the K~, K°,D~, D° axid B~ , B° mesons. Aj and 5i are the amplitudes 
with the two pions in isospin "I" states and their corresponding phases, respectively. 

The amplitude Aq is induced by an interaction with isospin 1/2, whereas A2 is induced 
by I = 3/2 interaction. The AI =1/2 rule refers to the fact that the ratio A2/AQ is much 
smaller than one. Using the above information one can write this ratio in terms of branching 
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TABLE I: Branching ratios. D — > 7T7t and B — > 7T7t are from Ref.|5j and Ref.[6|, respectively. 





Data 


naive factorization 


QCD factorization 


BR(D+ -> 7T+7T ) 


(2.6 ± 0.7) x 1(T 3 


2.97 x 10~ 3 


2.88 x 10~ 3 


BR(£>° -> 7T+7T-) 


(1.38 ±0.05) x 10~ 3 


2.02 x 10~ 3 


2.08 x 10~ 3 


BR(L>° -> vr ^ ) 


(8.4 ±2.2) x 10~ 4 


5.51 x 10~ 6 


3.65 x 10~ 6 


BR(L>° -► tt+tT + vr vr ) 


(2.22 ±0.23) x 10~ 3 


2.02 x 10~ 3 


2.08 x 10~ 3 


Rd 


0.65 ±0.13 


0.79 


0.75 


BR(5+ -> tt+vt ) 


(5.5 ±0.6) x 10" 6 


6.59 x 10~ 6 


5.84 x 10~ 6 


BR(£>° -> vr+vr-) 


(4.6 ± 0.4) x 10~ 6 


8.23 x 10" 6 


8.38 x 10~ 6 


BR(£° -> vr ^ ) 


(1.51 ±0.6) x 10~ 6 


2.29 x 10~ 7 


2.48 x 10~ 7 


BR(B° -> vr+vr" + 7r°vr ) 


(6.11 ±0.52) x 10~ 6 


8.46 x 10" 6 


8.63 x 10~ 6 




1.11 ±0.20 


0.96 


0.84 



The branching ratios for D — > nn have been measured to good precision^]. Recently 
the B — > 7T7T branching ratios have also been measured In Table IJ we list the branching 
ratios of the 7T7T decays of D and B. With these measured branching ratios, one can easily 
check if there are A J =1/2 rule in B and D to tttt decays. We obtain 



R D = 0.67 ±0.13; i? B = 1.11 ± 0.20. (3) 

It is obvious that the AI =1/2 rule is violated in B(D) — > tttt decays. The situation is very 
different from that in K — > tttt decays. It is therefore important if one can understand the 
situation by studying these decays more carefully. 

To have a full understanding of the situation, one also needs to consider phase shifts. 
From eq. (JTJ) one obtains the isospin phase differences cos 5p = cos(<5 — <5 2 ) p as 



3|A+-| 2 -6|yl 00 | 2 ±2U- | 2 
cos 5 P = — = — . (4) 



4V3\A-°\ ^\A+~\ 2 + |A 00 | 2 - 2|A-°| 2 /3' 



3 



Using experimental data, we obtain 

cos 5 D = 0.13 ±0.06, for D — ► 7r7r, 

cos5 B = 0.58 ± 0.20, for B -> tttt. (5) 

These results indicate that the difference of the phase shifts 8p is sizeable. In the above we 
have neglected CP violation in the decay amplitude which is applicable for D — > titt decays. 
If CP violation is sizeable which may happen in B — > decays, one needs to replace cos<5p 
by cos(5p + <pp) with cf>p being the CP violating phase difference in the decay amplitudes. 
We will come back to this later. 

We first discuss the quantity Rp. In the past few years progresses have been made in the 
calculations of a heavy meson decays into two light mesons based on QCD factorization. Sev- 
eral ways of calculating B — > inr have been developed with different methods treating non- 
purterbative quantities involved 0,0,0,13]. There are also many model-independent studies 
about the amplitudes for B — > 7T7T decays based on isospin and/or SU(3) symmetries jl3. 14 1. 
We will take QCD improved factorization as the theory to compare with data. We obtain 
the isospin amplitudes in the following 



r) In p p 
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-\/^'ubi + (Al + A c )(fo 3 + 26 4 - + -V)]^> 
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A 2 = - — [\' u (a 1 + a 2 ) + -X' p (-a 7 + -fxa p 8 + a 9 + a p w )}A 7Tn . (6) 

In the above expressions, = V qb V* d , A n7T = i(G F /\^2)(rn 2 B — ml)F^' 7T (0)f n , and Jx(^) = 
2m^/mi,(/i)(m M (/i) +m l i(fi)). The other quantities are defined in Ref.[3] 

Terms proportional to bi in eq. (JHJ) are referred as annihilation contributions. These terms 
have end-point divergences of the form Xa = f$ (j) a (y)dy/(l — y) and need regularization. 
The term proportional to f 11 also has a similar end-point divergence X^- These divergences 
associated with 1/mhn corrections indicate in a way the incompleteness of factorization 
calculation. In Ref.|7( these end-point divergences are parameterized as \ia(m b / A h )(l + 
Pa,h^ a,h )- An idea of the size of the corrections can be obtained by varying pa,h and 

<j>A,H- 

To finally obtain numerical results, one needs to know the CKM and the hadronic param- 
eters in the above amplitudes. There are considerable progresses in the determination of the 
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CKM parameters 0, Q|. We will use the central values given in Ref.jsJ with s±2 = 0.2243, 
S23 = 0.0413, Sis = 0.0037 and the CP violating phase 7(^13) = 60° for illustration. For 
the hadronic parameters we use the "default" values given in Ref.|?J- The results for the 
ratio Rb are listed in Table HJ In the table, we also list the values obtained by using naive 
factorization [12| for comparison. We see that the a s order correction is at the order of 10% 
and can be as large as 30% if the involved hadronic parameters vary within a reasonable 
range. One also notes that the term // generates an absorptive part which is absent in naive 
factorization calculations. 

One can see from Table |J that the theoretical calculation is in qualitative agreement with 
the data that there is no A J = 1/2 rule in B — > tttt decays. This is expected since that 
in factorization calculations the isospin does not play a special role. The leading operators, 
i.e. the tree operators, which contribute to B — > tttt, contain both AI = 1/2 and 3/2 pieces 
with similar weights, the amplitudes for isospin I =0 and I =2 are therefore expected to 
have similar sizes. 

We now consider the situation for D — > tttt decays. In Ref. Sannino noticed the 
violation of the AI = 1/2 rule in D — > tttt and tried to understand it based on the effective 
weak Hamiltonian. Assuming that the hadronic matrix elements M 2 , M Q and M defined in 
Ref. [l5 1 are all equal, numerically Rd ~ 0.32 ~ 0.44 which is not too far away from data. 

The c-quark is not as heavy as the b-quark, factorization calculation for D — > tttt may 
not work as well as for B — > tttt. We however expect that a theoretical calculation based 
on QCD factorization scheme can still provide some crude estimate. We therefore also use 
the QCD factorization for D — > tttt. In this calculation we will use the Wilson coefficients 
given in Ref.j3]. We have also checked the dependence on the renormalization scale and 
find that the result on Rd does not change much with respect to the scales. The numerical 
results are listed in Table |1J We see that the situation is similar to that for B — > tttt, there 
is no A J = 1/2 rule. 

The calculation for K — > tttt using factorization becomes very questionable. Nevertheless, 
attempts have been made to evaluate the amplitudes. Buras et al. evaluated the ratio 
Rk based on the large N expansion approach with 

R 1 A ci{p) < giQf) >o + c 2p) < 02<X) >o +Ei=3 C »(/ i ) < Oiif*) >o 
K A 2 ci(/i) < Oi(/x) > 2 +c 2 (/i) < 2 (//) > 2 
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and the hadronic matrix elements can be parameterized as 

<O 1>0 =- l -XB^\ <O 2 > = 5 -XB^\ 

< 3 > = \xB^ 2 \ < 4 > =< 3 > + < 2 > - < 1 > , 

< 5 > = \ < 6 > , < 1 > 2 =< 2 > 2 = ^XB[ 3/2) , 

< O e >„= -iJl[ _ 1^4"", (8) 

V 2 m s (/x) + m^/i) k 

where X = \J\F^{m 2 K — m%), k = F w / (F K — F n ) and Bi are the hadronic parameters to be 
determined either by fitting experimental data or invoking phenomenological models. By 

(3/2) (1/2) 

fitting data they obtained a set of the hadronic parameters as B\ 1 ' — 0.48, B\ = 10, 
B^ 2 ^ = 5 and B^ 2 ^ = B^ 2 ^ = 1. The value of Rk calculated in this approach agrees 
with data. One notes that several bag parameters Bi are substantially away from the naive 
value of Bi = 1. Thus one expects that some non-perturbative effects and the final state 
interaction effects are involved altogether in the parameters. 

We note that although theoretical calculations agree with data, one can be confirmed 
with the fact that there is not a AI = 1/2 rule in B(D) — > nn decays, there is a large 
difference in B°(D°) — > 7r°7r° branching ratio between naive theoretical estimates and data. 
It is important to see if these discrepancies can be explained. 

A possible source may be due to uncertainties in the factorization calculation, in particular 
corrections of l/m,b tC . As pointed out earlier that the end-point divergences appearing in the 
hard scattering and annihilation signal incompleteness of the l/m b c corrections. A complete 
treatment of these effects are beyond the scope of this work. To have some idea about the 
effects of l/m^ c we have carried out a calculation by varying the parameters pah from to 3 
and 4>a,h from to 2n. We find that the changes on the branching ratios of D + {B + ) — > 7r + 7r° 
and D°(B°) -> tt+tt- are in the range of 20% to 30%, and the changes on D°(B°) -> 7r°7r° 
can be dramatic (a factor of 2 to 3). However, it is still not possible to bring D°(B°) — > 7r°7r° 
branching ratios to their experimental values. One may need to consider other effects which 
are not included in the factorization approximation. 

To this end we consider long distant final state interaction effects. If the differences are 
really due to FSI effects, it should not only explain the differences mentioned above, but 
should also explain why there is a AI = 1/2 rule in K — > nn, but not in B(D) — > nn. 
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To explain the AI = 1/2 rule in K-decays, Isgur et al. [2J proposed a possible mechanism 
for the smallness of Rk (~ 1/22) that before the two pions produced from the weak decay 
of the kaon fly apart to become free particles, they reside in bound states of either I = or 
1 = 2. Due to different interactions for the two different isospin states, the wavefunctions 
at origin for 1 = and 1 = 2 states are distorted differently. The hadronic matrix elements 
of the two isospin states would undergo an enhancement or a suppression as [2] 

< (7T7r) J=0 I H^ 2 | K > dp i < (tttQS I H]l 2 \ K > 

< (vnr) /=2 | E% 2 | K > V < (tttOS | H% 2 \K> 

Here the distorting parameter di is defined as dj = \ip true (0)/ip^ ree (0)\ 2 . 

The general requirement for the interaction potential is that the 1 = channel experiences 
an attractive interaction whereas the 1 = 2 channel experiences a repulsive one. We note that 
a simple square well potential of the form XJ\\ ■ 1 2 = (1/2) U[I(I + 1) — Ji(Ji + 1) — l2(h + 1)] 
with U > can result in the required potential form. 

Fitting data on the low-energy 71 — 71 scattering phase shifts, the parameters of Vj and 
aj can be determined. The distortion due the potential can result in an enhancement of the 
channel of the AI = 1/2 over the AI = 3/2 by a factor r = (d /d 2 ) 1/2 for K — > 7T7T which 
can be as large as 9 ~ 10. For the isospin correlated potential UI\ ■ I2 with U = 0.4 GeV 
and potential range ao = a 2 = 0.8 fm, the enhancement factor r is 9. This factor can make 
up the gap between theory and experimental data leading to the AI = 1/2 rule in the case 
of K — > 7T7T. One can also try more complicated potential forms, such as Gaussian- type 
potentials as studied in Ref.j2|, but the qualitative features are unchanged. 

Since the final two pions produced from B and D decays have similar isospin structure 
except that their kinetic energies are much higher, they should experience similar effects 
due to the potential. However, since the two pions have much higher kinetic energies they 
have much shorter time to interact with each other before emerging out of the bound states 
to produce large effects. More specifically, in the D and B cases, the momenta |k| = 
\\J {M"^^ — M%) are much larger than the potential energies \Vi\, one can easily verify 
that do w d% ~ 1. Even though for the B and D decays, the two pions fly very fast and 
the non-relativistic quantum mechanical scenario may not produce accurate numbers, the 
qualitative conclusion about smallness of the bound state effects, however, should hold. 
We therefore see that this mechanism can easily explain why there is a AI = 1/2 rule for 
K — > 7T7T but not for B(D) — > tttt. 



We now discuss the FSI phases. Using experimental data, we obtain cos<5d and cos<5d 
which are shown in eq.(5). These results indicate that the difference of the phase shifts Sp 
is sizeable. Although factorization amplitudes in eq.© have included some re-scattering 
effects at quark level, they are not large enough to explain data. Factorization calculations 
using the default values for relevant hadronic parameters would give a cos Sp to be very 
close to 1. There are uncertainties in factorization calculations in regularizing the end-point 
divergences. As mentioned earlier, we have checked numerically that within reasonably 
allowed parameter space, these contributions cannot generate large enough phase shifts to 
reproduce data. 

The distorting effects in eq.fjHJ) discussed earlier can also generate a phase shift. At the 
kaon energy, it is possible to generate the required phase shift for K — > tttt. However these 
phase shifts go to zero as the energies become much larger than the potentials. For B and D 
decays, the phase shifts are practically zero because the pions from D(B) decays have high 
energy and pass through the potential in too short a time to produce any significant phase 
shifts. There is a need of additional effects. The effects producing large phase shifts in D 
and B decays must be different from the mechanism in Ref.Q which is indeed available, so 
that at higher energies more channels become active in producing absorptive part of decay 
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amplitudes and therefore phase shifts 

From the calculated branching ratios of D° — > tttt and B° — > tttt in Table HI we note that 
the value for tt + tt~ is obviously higher than data, whereas the value for 7r°7r° is significantly 
lower than data. However the sum of BR(P° — > tt + tt~) and BR(P° — > 7r°7r°) , and also 
BR(P~ — > 7T~7r°) are close to experimental values in both D and B cases. Since only these 
sums determine the ratios Rp (see eq.(2)), the theoretical calculations are close to data. One 
should keep in mind that there are uncertainties in many of the hadronic parameters. The 
theoretical predictions for the branching ratios can change, but it is difficult to generate 
large enough branching ratios for B°(D°) — > 7r°7r°. If re-scattering of the two pions can 
change the higher valued tt + tt~ into the lower valued 7r°7r° after the pions come out of the 
bound states as discussed above, it may be able to produce the correct values to meet data. 
In fact large phase shifts due to re-scattering in B(D) decay into two light mesons have 
been noticed beforell3l Il4. We have carried out a simple exercise similar to some work 

n 

in Ref.[13] taking the magnitudes of A 2 as determined by factorization calculations and 
S as a free parameter to fit data. We find that with Sd equal to 0.457T, cos<5d is close 



8 



to the central value of the data. In the case of B — > 7r7r, with the default values for the 
hadronic parameters, 5b = 0.307T can reproduce a cos 5b which makes the branching ratio 
of B — > 7r°7r° to be close to the experimental central value. In this case both amplitudes 
A 0t 2 are slightly higher than data, but can easily be made to agree with data if a slightly 
smaller form factor F B ^ n = 0.23 is used instead of the default value of 0.28. A more precise 
theoretical evaluation of the re-scattering phases would be difficult because the mechanisms 
are not fully understood yet. 

We finally make some comments about isospin analysis and CP violation. The FSI 
phases are very important for the study of CP violations. Sizeable CP asymmetry can occur 
in B — > Tin decays. We therefore will concentrate on B — > txtx decays. The isospin amplitude 
given before can be further decomposed into different components according to the CKM 
matrix elements associated with 

Ai = V ub V: d aJ + V tb V t * d a p . (10) 

The components aj mainly come from tree amplitudes by exchange W-boson but also small 
corrections from loop-induced c and u penguin contributions, and af are induced at loop 
level by c and t penguins. 

In general there can be 4 complex hadronic parameters 2 and a p 2 to describe B — > ixix 
decays. Among them one can always set one of the components to be real, there are actually 
only 7 independent parameters. Furthermore in the SM, the leading order from c and t 
penguin contributions to of ate dominated by t penguin front opetatot O 9 , 10 which have the 
same Lorentz structure as the tree operators. One has(lS| a P = 3/2(cg + C\o)/(c\ + C2)a T . 
Neglecting other smaller contributions to a 2 ' components, one only needs five hadronic 
parameters to describe B —>■ ixix decays. This is an important fact which has interesting 
implications. 

One of them is that it can be used to determine CP violating phase 7 if the 7 possible 
experimental observables in B — > tttt decays are all measured (the 3 branching ratios, the 
2 direct CP violation (B~ — > 7t~-k° has very small CP violation in the SM and probably 
cannot be measured), and 2 mixing induced CP violation for B° — > 7r + 7r _ , 7r~7r ), the 7 
observables can determine the five hadronic parameters and two CKM parameters p and 77. 
We need to wait for more data to carry out a full analysis. If one takes the CKM parameters 
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as determined by other data, the hadronic parameters can already be extracted. 

In obtaining the value cos Sp in eq.(jHJ, we have neglected possible CP violating effect in 
the isospin amplitudes. We therefore cannot determine and a<^ separately. One needs 
more data points to obtain more detailed information about the size of the decay amplitudes 
and also phases. Fortunately experimentally, besides the branching ratios listed in Table 
for B — > 7T7T, there are also some measurements on the direct CP violating parameter Aqp 
and mixing induced CP violating parameter Sep fj, 

A C p(B° -> 7T+7T-) = 0.37 ± 0.11, S CP (B° -> 7T + 7T~) = -0.61 ± 0.14 

A CP (B° -> 7r°7r°) = 0.28 ± 0.39. (11) 

Using five data points, the three branching ratios and the CP asymmetry parameters 
Acp and Sep for B° — > 7r + 7r~, as input we can solve for the hadronic amplitudes. We have 



1) 4 = 0.366,4, = 0.276e- i48 ' 5 °A a£ = 0.0Q7e~ i71AO A; 

2) 4 = 0.366A, = 0.289e i54 - 3 °A, a£ = 0.065e- il66 ° A (12) 

Here A = ^BR{B~ — > 7r vr )rs toia , Factorization calculations give different values, in par- 
ticular the phases. Improved theoretical framework of calculating the amplitudes is needed. 

If the penguin amplitude Aq is neglected, there is no direct CP violation, one can use eq. 
P| to determine the FSI phase. Since Aq is small compared with the three amplitudes, the 
use of eq. (0) is a good approximation. We note that the two solutions have different FSI 
phases for each individual amplitude, but the phase difference of Aq and A\ are similar in 
size for solutions 1) and 2) and different in signs, Since the cosine function is not sensitive 
to the sign, one would obtain similar cos<5b. One needs to find some observables which can 
distinguish these solutions. We indeed find that CP violation in Acp of B° — > 7r°7r° depends 
on the sign. We have not used the data as input because the error is large there. Using 
the above amplitudes and phase we can predict its value. We obtain Acpij ^ ) = —0.60 
and 0.18 for solutions 1) and 2), respectively. The solution 2) obtains a value close to the 
current central value. But due to the large error at present, it is too early to decide which 
solution is the correct one. When more precise data become available one can distinguish 
these solutions and obtain more detailed information about the isospin amplitudes. 
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In summary in this work we have analyzed the isospin amplitudes due to AI = 1/2 and 
AI = 3/2 interactions in the K — > tttt, D — > tttt and B — > 7r7r decays. Experimental data 
clearly show that unlike the situation for K — ■> 7T7T, there is no AI = 1/2 rule in B(D) — > 7r7r 
decays. Theoretical calculations using the naive factorization and the QCD factorization 
are consistent with data about the violation of AI = 1/2 rule, but there are difficulties to 
obtain correct phases in the amplitudes and branching ratio for BR(B° — > 7r°7r°). 

The question why there is a A/ =1/2 rule in K — > 7T7t whereas not for B(D) — > 7T7t 
can be interpreted by the re-scattering of the final two pions due to a isospin correlated 
interaction proposed by Isgur et al. We note that the potential of the form {71 1 • I 2 = 
(l/2){7[/(7 + 1) - Ji(/i + 1) - J 2 (-?2 + 1)] can Piay such a role. To fully explain B(D) -> tttt 
data, there is also the need of re-scattering after the pions come out of the bound states. 
These additional phases are not precisely calculable at present. 

We have also shown that isospin analysis for B — > tttt provides valuable information about 
CP violation in B decays and hadronic amplitudes. 
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